function [ LL ] = PoissLambdaLogL(Lambda,Data)
%PoissLambdaLogL Log Likelihood of lambda given Poisson distribution
%   Lambda:     Expected counts
%   Data:       Set of observed counts

%This is the log-likehood of the Poisson distribution using 
%The Stirling approximation for log(k!)
LL=sum( Lambda-Data.*log(Lambda)+(Data.*log(Data)-Data));


end

